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ABSTRACT 

Context. The orbital inclination of the symbiotic prototype Z And has not been established yet. At present, two very different values 
are considered, ; ~ 44° and ; ^ 73°. The correct value of i is a key parameter in, for example, modeling the highly-collimated jets of 
Z And. 

Aims. To measure the orbital inclination of Z And. 

Methods. First, we derive the hydrogen column density (/jh). which causes the Rayleigh scattering of the far-UV spectrum at the 
orbital phase (f> = 0.961 ± 0.018. Second, we calculate ;ih as a function of ; and for the ionization structure during the quiescent 
phase. Third, we compare the nnC^ 0) models with the observed value. 

Results. The most probable shaping of the H i/H ii boundaries and the uncertainties in the orbital phase limit ; of Z And to 59° - 
2° / + 3°. Systematic errors given by using different wind velocity laws can increase / up to ~ 74°. A high value of ; is supported 
independently by the orbitally related variation in the far-UV continuum and the obscuration of the Oi] /11641 A emission line around 
the inferior conjunction of the giant. 

Conclusions. The derived value of the inclination of the Z And orbital plane allows treating satellite components of Ha and H/3 
emission lines as highly-collimated jets. 

Key words, binaries: symbiotic - Scattering - Stars: individual: Z And 



1. Introduction 

Symbiotic stars are long-period interacting binaries comprising 
a cool giant as the donor star and a compact star, most often a 
white dwarf (WD), as the accretor This composition requires 
large orbital periods, which are typically of a few years. 

The circumstellar environment of symbiotic stars, which 
consists of energetic photons from the hot star and neutral par- 
ticles from the cool giant, represents an ideal medium for the 
effects of R ayleigh and Raman scatte ring on neutral hydrogen 
atoms (e.g.. Nussbau mer et al. II1989I) . The Rayleigh scattering 
produces a strong attenuation of the continuum around the Ly-a 
line in systems with a high orbital inclinati on at positions wit h 
the giant in front during quiescent phases dlsliker et al. II1989I) . 
The strength of this attenuation is thus determined by mh be- 
tween the hot star and the observer Its quantity is a function of 
the binary position and the orbital inclination, depending on the 
ionization structure of the binary. For the steady state situation 
and the q uiescent phase, the Hi/Hii boundaries were first de- 
scribed bv lSeaquist et al l ([1984) (hereafter STB). In the model, 
the neutral zone is located symmetrically around the binary axis 
and is usually cone-shaped w ith the giant at its top facing the 
ionizing source (e.g.. Fig. 7 of lFernandez-Castro et al. 1 [T988l) . 

Z And is the prototype of the symbiotic stars class. The 
binary consists of a late-type, M4.5III giant and a WD ac- 
creting from the giant's wind on the 759-da y orbit (e.g., 
iNussbaumer & Vogel I [19891 iFekel et alTI l2000bl) . Z And was 
conside red to be a non-echpsing binary. B ased on optical obser- 
vations, iMikolajewska & Kenvon I ([1996) suggested the orbital 
inclination i » 50°-70°. Based o n polarimetric mea s uremen ts of 
the Raman scattere d /16830 line. ISchiiild &SchiIdl(fl997al) and 
llsogai et al. I(l2010l) determined that ; = 47°+12° and41°+8°,re- 



spectively. On the other hand, a high orbital inclination was sup- 
ported by (i) the eclipse-like effect measured in the light curve 
at the position of th e inferior conj unction of the giant, which 
suggested / > 76° jSkopal I l2003h . (ii) the polarimetric mea- 
surements in the continuum, whose modul ation along the orbit 
con-esponds to / = 73° ± 14° (llsogai et al. '"2010). (iii) the two- 
temperature type of the ultraviolet (UV) spectrum that develops 
during active p hases of systems with a high orbital inclination 
(ISkoDalll2005l i. and (iv) the presence of the Rayleigh attenuated 
far-UV continuum measured around the inferior conju nction of 
the giant during the quiescent phase (e.g. lSkopal II2003I) . 



Knowledge of the correct value of ; is a key parameter 
in, for example, modeling th e enhanced bipolar mass outflow, 
launched recendy by Z And dSkopal & Pribullall2006l) . If ; < 
55°, t his outflow c an be treated as a radiat ively accelerated 
wind dTomov et al.|[2010 : Kil pio et all 1201 lb . Otherwise, this 
type of mass outfl ow represents highly colhmated jets (e.g., 
ISkoDaIetal.ll2009h . 



In this paper, we model the effect of the Rayleigh scattering 
measured in the spectrum of Z And from the quiescent phase 
to derive /. In Sect. 2, we introduce observations and derive the 
riu parameter, which causes the attenuation effect. In Sect. 3, we 
introduce the ionization structure of the binary and calculate «h 
as a function of the orbital phase and the inclination. In Sect. 4, 
we compare the nu{i, <f>) models with the measured value and 
in this way derive the range of / from Rayleigh scattering. The 
conclusions are found in Sect. 5. 
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Fig. 1. UV observations of Z And carried out close to the in- 
ferior conjunction of the giant during its quiescent phase (on 
February 3rd 1988). Compared is the model SED (heavy solid 
line) with its components from the nebula (dashed line) and the 
hot star (thin solid line). The strong attenuation of the contin- 
uum around the Ly-a line is caused by the Rayleigh scattering 
process. Dotted line represents the non-attenuated light (see text 
for details). 
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Fig. 2. A scheme of a symbiotic binary to calculate hh along 
the line of sight, /, throughout the neutral Hi region (Eq. (|4]i). 
The thick solid line represents the Hi/Hii boundary (Sect. 3.3) 
and 1? is the angle between I and the binary axis. Other auxiliary 
parameters are introduced in Sect. 3.1. Distances are in units of 
the binary components separation, p. 



2. Observations 

There are two pairs of UV spectra of Z And in the International 
Ultraviolet Explorer (lUE) archive taken close to the position 
of the inferior conjunction of the giant during the quiescent 
phase: SWP07410 + LWR06393, from December 15th 1979, 
and SWP32845 + LWP12624, from February 3rd 1988. Both 
spectra display signatures of the Rayleigh scattering attenuation 
in the continuum around the Ly-o' line. However, the former has 
underexposed the far-UV continuum and is affected by a strong 
geocoronal Ly-ff emission. Therefore, we treated only the lat- 
ter s pectrum. According to the orbital elements of iFekel et al. I 
(|2000b), the spectrum was obtained at the orbital phase - 
0.961 ± 0.018. Its spectral energy distribution (SED ) in t he con - 
tinuum was modeled in the same way as by Skopal I (l2005h . 
Here, we simplified the model by using only one nebula and 
adopting zero He^^ abundance (Fig.lTJ. The strong depression of 
the continuum around Ly-a is caused by the Rayleigh scattering 
of the continuum photons on the neutral atoms of hy drogen (e.g., 
Ilsliker et al. ll989HVogelll991l;lDumm et al. 119991) . In our case, 
the depression corresponds to n^ = (3.9 + 0.5) x 10^^ cm"^. 

Having the quantity of ne from observations at the given or- 
bital phase, we can calculate its dependence on the orbital incli- 
nation, i, for the ionization structure during the quiescent phase. 
We introduce the model in the following section. 



3. The model 

Here, we determine the ionization structure of a symbiotic binary 
and calculate the column density of neutral hydrogen on the line 
of sight as a function of (p and /. Comparing theoretical values 
of nii(i, (p) with those derived independently from observations 
allows us to estimate the inclination of the orbital plane. 

3.1. hh as a function of the orbital inclination 

Assuming that the wind from the red giant is spherically sym- 
metric and neutral, the continuity equation determines the parti- 



cle concentration of hydrogen as 



M 



47Tr^lJt7iiiv(r)' 



(1) 



where M is the mass-loss rate from the giant, r is the distance 
from its center, p. the mean molecular weight, mu the mass of the 
hydrogen atom, and v{r) velocity of the wind. Accordingly, we 
obtain the hydrogen column density, «h, by integrating Eq. ^ 
along the line of sight, I, from the observer (-oo) to the ionization 
boundary as 



M 



nn = 



4npm 



H J" 



dl 



(2) 



where is a segment of the line of sight calculated from the in- 
tersection of its normal, which connects the giant's center, to the 
boundary (Fig.|2]i. To determine nu we need to know the radius- 
vector r at each point of the line of sight, which is a function of / 
and the phase angle, (p{(p = Oat the giant's inferior conjunction). 
Furthermore, we introduce a parameter b as 



b = p (cos / + sin" tp sin 0. 



(3) 



which represents projection of the separation p between the bi- 
nary components into the plane perpendicular to the line of sight 
(Fig.|2]i. Equation (|2]) may now be expressed as 



M 



nn = 



r dl 

J (P + b^)v( \ 



AnpniH J (p + h2-^v( V/2 + b^) 



(4) 



Then, knowing the value of l^ for a given direction to the ioniz- 
ing source and the velocity v(r), we can solve equation (|4]i for «h- 
Figure |2] shows that l^ - ^p'^ — b^ - s^, where is the distance 
from the hot star to the H I/H II boundary on the line of sight. 
Its value is given by the ionization structure during the quiescent 
phase (Sect. 3.3), which requires knowledge of the velocity v(r) 
of the giant's wind, which we introduce below. 
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Fig. 3. H i/H n boundaries in a symbiotic binary given by solu- 
tion of Eq. ^ for the velocity law d?) of the giant wind. This 
represents basic ionization structure of the circumstellar matter 
in symbiotic binaries during quiescent phases as originally out- 
lined by STB (Sect. 1). 



3.2. Velocity profile of the wind from the giant 



IVogell ([r991 ) recognized that the column densities riu that he 
measured from the spectra of EG And cannot be explained by a 
standard ^S-law wind (see, e.g., his Eq. (6)). The observed rapid 
increase of «h values prior to the total eclipse requires a lower 
velocity in the vicinity of the cool giant (r ^ Rg) and a much 
steeper velocity profile for 0.09 <: <p ^ 0.14 {r ^ 3Rg with re- 
spect to that given by the ,6-law wind (see Figs. 4 and 7 of lVogel I 
1199 lb . Therefore. PVogel I (1199 lb employed Eq. (4) to reconstruct 
a more appropriate velocity field of the wind from cool giants. 

Having many more observations coveri ng the echpse profile 
of SY Mus due to the R ayleigh scattering. iPereira et al. I j 19991) 
and'Dumm e t al. I (11999) solved Eq. (4) for the velocity law v(r) 
by expanding the function nuib) into a Taylor series. To match 
the «H values measured around the eclipse, both groups of au- 
thors confirmed independently the claim of a significantly faster 
acceleration of the wind particles from the giant in comparison 
with the y6-law. 

Re-analyzing the wind acceleration zone for the giant in 
EG And by using a se ries of FUSE and HST/STIS observations, 
ICrowlev et al. I (l2005h recently confirmed that the wind acceler- 
ation region extends only to ~ 2.5 R„ from the limb of the giant. 
The correspondin g wind velocity profile was very similar to that 
derived bv ,Dumm et al. I ( IT999,) for SY Mus (see their Fig. 7). 
Therefore, for the purp ose of this paper, w e use the wind ve- 
locity law as derived bv lDumm et al. 1 (fT999h . We note that cool 
giants in both SY Mus an d Z And have the same spectral type 
(M4.5, Miirs et & Schm id 1999) an d sim ilar other fundamental 
parameters (see Table 2 of Skopal Il2005h . which also supports 
adoption of the w ind law (O for the ca se of Z And. 

According to ' Dumm et al. I (11999 *). it was sufficient to ap- 
proximate the nu(b) function by two terms of the Taylor expan- 
sion (k = I and k ^ 20), which allows expressing the velocity 
law in a form 



rv(r) Air Aicr''' 



(5) 



where a - 2M lAn^niHVcxRg, ai, Ok and k are parameters, and 
is defined recursively as 



7r/2, 



Ak^K/[(2k-2)At-i]. 



(6) 



Fitti ng the obse r ved n -g as a function of the impact parameter 
b. .Dummet all (1 19991) derived ai ^ 4.5 x lO^^cm a ^ 3 x 
10^^ cm~^, ^ =s 20 and a2o ~ 1 x 10-" cm"^. Using these values 
we can rewrite Eq. (|5]l in the form 



v(r) = 



1 + 



lO'* 



(7) 



3A2or 



19 



where r is the distance from the center of the giant in units of 
its radius. According to this law, the main acceleration zone of 
the wind is located between the giant's photosphere and r = 3 
(y(3)/yoo ~ 0.9, A20 = 0.2838). 

3.3. Ionization structure 

The distance s^p from the ionization boundary to the WD (Fig. [2]) 
is determined by a balance between the ionizing photons emit- 
ted in a small angle arou nd the direction -& and recomb inations in 
this angle. According to lNussbaumer & Vogell(ll987h . the equi- 
librium condition can be expressed as 



L^h— = — J NHAs)NAs)aB(H,TM7^s^ds, (8) 



where Lph is the number of photons capable of ionizing hydro- 
gen that are emitted spherically-symmetrically from the hot star 
per second, s is the distance from the hot star, Nh^ and A^^, are 
concentrations of protons and electrons, respectively, and ob is 
the total hydrogen recombination coefficient for recombinations 
other than to the ground state. In the region, we consider that 
Ne = Nh = Nh^ . Then, assuming that T^, is constant throughout 
the nebula and thus also a^, condition ([8]) can be expressed as 



where X^^ is the ionization parameter (STB), 



ffB(H, Te) 



and 



f — 



(r) 



ds. 



(9) 



(10) 



(11) 



where we used the particle concentration in the giant's wind as 
given by Eq. ([T]i, and 



- + - 2sp cos 
Examples of H i/H 11 boundaries are depicted in Fig. |3] 



(12) 
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4. Orbital inclination of Z And 

Having defined ionization structure, we can solve Eq. (|9]l for s^, 
which determines and thus hh (see Fig. |2] Eq. (4)) for the 
ionization parameter X^^ . Its value is given by the binary pa- 
rameters. Therefore, first we find a possible rang e of X^^ . 

For foib = 759 days ('Fekel et al. ' '2000H) and th e tota l 
mass of the binary, 2.6 Mq (Mikolaiewska & Ken vonI I l996h . 
we get the separation of the stars, p - A%QRq. Furthermore, 
we adopt aR(H, T,) ^ 1-43 x 10^'^cm^s"' for = 20000/: 
jNussbaumer & Vogell [T987h and M ~ 7 x lO'^Moyr"', 
deriv ed from the total emission measure during quiescence 
(ISko pal 2005). This value of M represents an upper limit be- 
cause of the presence of other sources of the nebular radia- 
tion, e.g., the hot star wind. The hot star luminosity Lj, ~ 
2 3 00tf/1.5kpc)^L fT-, and a temperature of -120 000 K (Table 3 
of ISkopall I200I yield Lph ~ 1.5 x lO^V/l-S kpc)^ s '. 
Assuming a typical val ue for the termin al velocity for the giant 
wind, ~ 40km s"' (Reimers 198_1) , we o btain X"^ x 20. 
We note that iFernandez-Castro et al. 



(Il988h derived X"^ 



14. However, there are more sources of uncertainties in deter- 
mining the parameter X'^^: (i) Terminal velocity for the mas- 
sive slow wind from the jiant was jilso assumed to be of 
20kms ' only (e.g.. lDumm et al. II1999I) . This value scales the 
above derived X^'^ with a factor of 0.25. (ii) A smaller dis- 
tance to Z And of 1.12kpc a nd 1 .2 kpc ( iFernandez-Castro et al. I 
ll988HSokoloski et al. l2006h reduces the luminosity and thus the 
value of X^"^ with a factor of 0.55 and 0.64, res pectively, (iii) 
A m aximum Zanstra temperature of 130 000 K jMiirset et al. I 
[199 1) gives Lh = 3 400(t// 1 .5kpc)2 L©, i.e. Lph = 2.2 X 
10'*^(d/1.5kpc)^s"', which enlarges the above value of X'^^ = 20 
with a factor of ~1.5. According to these possibilities we con- 
sider a range of the ionization parameter as 



3 <X"^ ^ 30. 



(13) 



In the following subsection, we explore the riuii, (p) function at 
the position of (p - 0.961 + 0.018, when the lUE observations 
were exposed. 

4.1. Modeling the column density at(p - 0.961 

The top panel of Fig.|4]shows the resulting column densities as a 
function of the orbital phase, calculated for the average value of 
X"+ = 17, for which nH(;, = 0.961 ±0.018) = 3.9 X 10^2 cm-2. 
The range of possible orbital phases limits the inclination an- 
gle to / = 59° - 2° / + 3°. In the bottom panel of Fig. H the 
nn(i, 4>) function was calculated for the range of possible X^^ 
(Eq. (fTTb). Here, the models for different values of X"^, satisfy- 
ing nii{i, 0.961) - 3.9 x 10^^ cm"^, restrict the orbital inclination 
to i = 59° -2° / + Small uncertainties reflect small differences 
between opening angles of H i zones for larger values of X"^ 
(see Fig. [3). Figure |5] illustrates the calculated nH(', 0.961) as a 
function of / for the range of the parameter X^^ . From the figure, 
we can see that the uncertainty in the measured value of hh can 
enlarge the resulting uncertainty in ; by a few degrees. However, 
systematic errors can result from using a different wind veloc- 
ity law. Although the used v(r) satisfies best the measured hh 
as a function of the parameter b (Sect. 3.2), its high values 
^ lO^'* cm"^ in the vicinity of the stell ar disk of the giant c annot 
be measured accurately (see Fig. 5 of iDumm et al. This 
precludes a correct determination of the wind law. The riuib) val- 
ues, as measured prior to and after the inferior conjunction of the 
giant, also reflect an asymmetry of the neutral wind zone w ith re- 
spect to the binary axis (see Fig. 4 of lDumm et 311119991) . The 
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Fig. 4. Top: Calculated column densities (Eq. (4)) as a function 
of the orbital phase for the H i zone defined by the mean value of 
- 17 (relations (10) and (13)). Particle concentration along 
the line of sight (Eq. (1)) was calculated for Vm = 40km s"'. 
Solutions giving «h = 3.9 x 10^~ cm"^ at the orbital phases = 
0.961 ± 0.018 (dotted lines) are plotted. Bottom: As in the top, 
but for 3<X"+ < 30 and S = 0.961. 



used wind law O was derived from the egress data. Thus, the 
steeper ingress data will correspond to a steeper v(r). Therefore, 
we also investigated a limiting case for a maximum steepness 
of v(r) = Voo. This approach yielded / = 74° .7 - 4°.0 / + 10°.3. 
From this point of view, the orbital inclination of Z And can be 
expected to be in the range of ~ 59° - 74°. 

Finally, the dependence of / on M follows that on X^^, 
because X"^ oc M^^ (Eq. (fTOt). This implies that M < up- 
per limit of 7 X lO"^M0yr"' (see above) yields higher X"^, 
which, however, has no significant effect on the corresponding 
value of /, because the opening angle of the neutral zone, i?oo, 
decreases very slowly for X^^ > 10 (see Fig. [3]l. For exam- 
ple, the lo wer value of M = 3.1 x 10"^(<i/1.5kpc)^''2 yr"', 
derived by iFernandez-Castro et al. ]([1988) from radio observa- 
tions, enlarges the range of X'^^ ( fT3T l by a factor of ~5, i.e. 
15 < X"^ ^ 150. Figure |5] also shows the model nH(;, 0.961) 
for X^^ - 100, which demonstrates only a small increase of / to 
60.°7. 
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Fig. 5. As in the bottom panel of Fig. IH but for hhC' , = 0.96 1 ) 
as a function of /. Solution for X^^ =100 was added to demon- 
strate the small change of / for very high X^"^ (Sect. 4.1). 
Intersections of the model curves with the horizontal line at 
the observed value of nu = 3.9 x 10^^ cm"^ correspond to 
; = 57, 59, 59.5, and 60. °7. The shadow band represents the 
uncertainty in the measured tin- 



4.2. The far-UV light curve 

The high orbital inclination of Z And is also supported by a large 
variation in the far-UV continuum along the orbit as measured 
on lUE spectra. Figure |6] shows examples of fluxes measured 
around /11280 A, which is relatively free of emission/absorption 
features. Similar variations are indicated also at longer wave- 
lengths, where the influence of the Rayleigh scattering is neg- 
ligible. For example, we selected fluxes at ~ /1 1600 A, where 
instead the iron curtain ab sorption could be pronounced (see 
IShore & AufdenberglfT993h . The variation here has a smaller 
amplitude with a shallow minimum with respect to that at ~ 
/11280A. Generally, the UV continuum peaks around the supe- 
rior conjunction of the giant, while at the opposite site the fluxes 
are significantly fainter. This variability reflects the presence of 
an additional source of extinction, which attenuates more the 
short wavelength part of the UV spectra. This effect is known 
for eclipsing symbiotic binaries RW Hya, S Y Mus, and EG And 
(iDumm etal. Ill 9991: ICrowlev et a"ni2008h The influence of the 



absorbing medium was modeled by Horne et alTI (Il994h . who 
found that the absorbing gas produces a modest optical depth 
in the Balmer and Pashen continuum, which can make it sig- 
nificantly lower in the UV (see their Fig. 8). As the absorbing 
material is concentrated more on the orbital plane and enhanced 
at positions with the giant in front, the large amplitude of the 
far-UV continuum can be attributed to a high orbital inclina- 
tion of Z And. In contrast, the non -eclipsing symbiotic binary 
AG Dra (; ^ 30 - 45° o r ; ~ 60°. iMikolaiewska et al. l[T99l 
ISchmid & SchildllI997bh does not show such variability in the 
UV spectrum (Fig.|6ll. 

Finally, we calculated the Rayleigh scattered light curves for 
the nii(i, (p) functions depicted in the bottom panel of Fig. 31 
scaled to the observed flux at = 0.961 (see Fig.|6l). The aim 
here is to model the expected Rayleigh attenuation as a func- 
tion of (f) and to check how it is connected with the neighboring 
measurements. Assuming that the Rayleigh scattering is the only 
process attenuating the far-UV continuum, we can write the ob- 
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Fig. 6. Variation in the far-UV, 1280 A (•) and 1600 A (+) 
fluxes of Z And and AG Dra along their orbits as measured 
on the lUE spectra during their quiescent phases. The arrow 
points the flux from the spectrum in Fig. [1] The lines around 
= and 1 represent Rayleigh attenuated light curves for mod- 
els in keys. Orbital elements of iFekel et al. i (.2000b.) were used. 



served flux as 



■0 „-o 



(14) 



where F° is the original flux and cr^'^^ is the Rayleigh scattering 
cross section. For Fi28o(0.961) = 6.5 x lO"'-' erg cm"^ s"' A"' 
(Fig. IS, «h(0-96 1) = 3.9 X 10^^ cm'^ (F ig. [D and crf^^^,, = 
1.1 X 10"^^ cm^ dNussbaumer et al. I Il989h we get the unscat- 
tered flux F^l2m - 1 -0 X 10^'^ erg cm~^ s"' A"' . Then, according 
to Eq. (fT4l i. we calculated the light curve Fi28o(0) for the three 
models nuii, (p) plotted in Fig. |4] The result is shown in Fig. |6| 
The rather narrow 0.1 Porb) and deep model light curves re- 
flect a close cut of the H i region for the resulting inclinations 
and the strong dependence of the Rayleigh attenuation on hh- In 
other words, due to a small opening angle of the neutral zone 
(i?oo ~ 35° and 34° for X = 17 and X = 30, respectively), 
their cut with the line connecting the hot star and the observer 
along the orbit covers only a small segment of the orbital period. 
Overall, the Rayleigh scattering light curves, calculated accord- 
ing to our resulting «H(i, (p) models, seem to be consistent with 
the fluxes measured around the inferior conjunction of the giant. 



4.3. A probe of a high inclination with the O i] 1641 A line 

IShore & Wahlgreni (|20IO|) investigated observational properties 
of the intercombination transition Oi] 1641 A to a metastable 
state in the spectra of some symbiotic stars with the aim of di- 
agnosing the neutral red giant wind. Their idea is based on the 
following: (i) the ionization potential of both oxygen and hydro- 
gen is comparable, (ii) the oxygen resonance multiplet (1302, 
1304, 1305 A) is, in addition to the ground state, connected to 
two metastabl e states through emissio n at 1641 A and 2324 A 
(see Fig. 1 of iShore & Wahlgrenll20IQi) . and (iii) the intercom- 



5 



A. Skopal and N. Shagatova: Measuring the orbital inclination of Z Andromedae from Rayleigh scattering 







1 i 1 1 1 1 1 1 

: ZAnd 

- i=60° n 

□□□ ° □ □ 

.B , , , □ , aim 


□t 


El 




□ 


□ 

□ □□ . 

B 


EG And □ 

- i=90° ^ 

□ □ ° □ □ 
1 . . . 




J 


□ 


□ 

□ 



0.5 1 

Orbital phase 



Fig. 7. Variations of the Oi] 1641 A line fluxes during quies- 
cent phases of Z And and EG And as a function of the or bital 
phase according to the ephemeris of iFekel et al. I (l2000al) and 
iFekel et aTl ( l2000b ). respectively. The light gray belts corre- 
spond to the intersection of the H i zone with the line of sight 
along the orbit (see Sect. 4.3). The dark belt reflects a time inter- 
val of the total eclipse of the point source by the giant in EG And 
(A0 = ±0.038; see Table 1 of Crowlev et al. 2008). Fluxes were 
taken from Figs. 4 and 8 of iShore & Wahlgreni (i2010i) . 

bination transition at /11641 A can be stimulated by absorption 
from the ground state of the 1 1302 A resonance line in the sur- 
rounding neutral gas. This implies that the Oi 1302 A emission 
line, which is formed by recombinations within the extended H ii 
region, surrounds closely the neutral H i region during quiescent 
phases of symbiotic stars and thus gives rise to the measurable 
O i] 1641 A emission within th e neutral wind from the cool gi- 
ant. IShore&Wahlgren I (120101) found that the O i] variations are 
strongly correlated with the optical light curves and also with the 
O VI Raman emission at ^16825, 7082 A. Both the O i] 1641 A and 
the Raman lines are optically thin; the transitions are not repeat- 
able: once they happen, their photons escape the medium. As 
a result, the transitions arise in the red giant wind close to the 
H i/H II interface, preferentially between the binary components 
with the highest density of the giant's wind (in contrast to the 
Rayleigh scattering, which is a repeatable pr ocess and thus mea- 
sures t he total thickness of the neutral wind). [Shore & Wahlgren I 
(l20Toh confirmed this view for the eclipsing system EG And, 
where the O i] line was obscured around the inferior conjunc- 
tions of the giant. For low-inclination systems, the authors sug- 
gested that these lines should always be visible. 

Accordingly, to probe the orbital inclination of Z And, we 
plotted the Oi] 1641 A line fluxes of Z And from its quies- 
cent phase together with those of the eclipsing symbiotic binary 
EG And in the phase diagram (Fig. |7]l. The figure demonstrates 
that around the inferior conjunction of the giant ((/> = or 1), the 
Oi] fluxes are very weak (EG And) or not measurable (Z And). 
Maximum values were measured when the line of sight does 
not intersect the neutral hydrogen zone (the light gray belts in 
Fig. |7]l. For EG And, the extent of the Hi region in the phase 
diagr am (A0 = +0.16) was determined for ; = 90° ( Vogel et al. 
1 19921) an d the position of the H i/H ii boundary at b 3.5 Rg (see 
Fig. 7 of ICrowlev et al. |[2005L and Eq. Q here). In the case of 



Z And, the range of intersections of the line of sight with the 
neuti-al zone in the phase diagram (A(p - ±0.11,0.048,0.044; 
dotted lines in Fig. |7]i corresponds to our models for X - 3, 17, 
30 from the bottom panel of Fig. |4] The results are consistent 
with the properties and the formation region of the Oi] 1641 A 
line during quiescent phases, as described above. Thus, the oc- 
cultation of the O i] fluxes of Z And during its quiescent phase 
independently supports a higher inclination of its orbit. 

5. Conclusion 

We used lUE observations of the symbiotic star Z And, which 
were taken during its quiescent phase around the inferior con- 
junction of the giant, to determine the inclination of its orbital 
plane. 

First, we derived ne - 3.9 + 0.5 x 10^^ cm"^ (Fig.[IJ on the 
line of sight, which passes throughout the neutral zone of the 
giant's wind towards the hot star at the time of the observation 
(at (p = 0.961 ± 0.018), by modeling the Rayleigh scattering of 
the far-UV continuum (Sect. 2). Second, we determined the STB 
ionization structure for the velo city profile of the giant's wind as 
derived by iDumm et al. ] (119991) (Fig. |3]l, and calculated nii{i, (p) 
as a function of the orbital inclination and the phase (Figs. |4] 
andjSj. Third, comparing the nii(i, (/>) function with the observed 
value of 3.9 x 10^^ cm"^, we estimated the range of orbital in- 
clinations as follows: (i) for the average value of the parameter 
X"^ - 11, i - 59° -2° /-i-3°, where the uncertainties reflect those 
in the binary position, and (ii) for the most probable position of 
the binary at = 0.961, the range of 3 :$ % 30 restricts 
the orbital inclination to / = 59° - 2° / + 1°. Systematic errors 
given by using different wind velocity laws can increase / up to 
~ 74° (Sect. 4.1). Fluxes, that are measured around the inferior 
conjunction of the giant are consistent with the Rayleigh scat- 
tering light curves (Eq. (fT4l l) calculated for our resulting iinii, (p) 
models (Fig.|6]l. 

The high value of the orbital inclination is also supported by 
the large amplitude of the orbitally related variation in the far- 
UV continuum (Sect. 4.2, Fig.|6ll and by the obscuration of the 
Oi] 1641 a fluxes around the inferior conjunctions of the giant 
during the quiescent phase of Z And (Sect. 4.3, Fig.|7]). 

Finally, the higher value of the inclination of the Z And 
orbital plane allows us to interpret the satellite components of 
Ha and H/3 emission lines, which appeared during active phases 
from 2006 as highly coUimated jets rather than as a special type 
of radiatively accelerated wind. 

Acknowledgements. We thank the anonymous referee for constructive com- 
ments. This research was supported by a grant of the Slovak Academy of 
Sciences, VEGA No. 2/0038/10, and by the realization of the Project ITMS No. 
26220120029, based on the supporting operational Research and development 
program financed from the European Regional Development Fund. 

References 

Crowley, C, Espey, B. R., & McCandhss, S. R. 2005, In Proc. of the 13th 
Cambridge Workshop on Cool Stars, Stellar Systems and the Sun, eds. F. 
Favata, G. Hussain & B. Battrick, (ESA SP-560: ESA), p. 343 
Crowley, C, Espey, B. R., McCandliss, S. R. 2008, ApJ, 675, 711 
Dumm T., Schmutz W., Schild H., Nussbaumer H., 1999, A&A, 349, 169 
Femandez-Castro, T., Cassatella, A., Gimenez, A. Viotti, R. 1988, ApJ, 324, 
1016 

Fekel, R C, Joyce, R. R., Hinkle, K. H., & Skrutskie, M. 2000a, AJ, 119, 1375 
Fekel, R C, Hinkle, K. H., Joyce, R. R., & Skrutskie, M. 2000b, AJ, 120, 3255 
Home, K., Marsh, T. R., Cheng, F. H., Hubeny, I., & Lanz, T. 1994, ApJ, 426, 
294 

Isliker, H., Nussbaumer, H., & Vogel, M. 1989, A&A, 219, 271 

Isogai, M., Seki, M., Ykeda, Y, Akitaya, H., Kawabata, K.S., 2010, AJ, 140, 235 



6 



A. Skopal and N. Shagatova: Measuring the orbital inclination of Z Andromedae from Rayleigh scattering 



Kilpio, E., Bisikalo, D., Tomov, N., & Tomova, M. 201 1, Ap&SS, 335, 155 
Mikolajewska, J., Kenyon, S. J., Mikolajewski, M., Garcia, M. R., & Polidan, R. 

S. 1995, AJ, 109, 1289 
Mikolajewska, J., Kenyon, S. J., 1996, AJ 112, 1659 

Murset, U., Nussbaumer, H., Schmid, H. M., & Vogel, M. 1991, A&A, 248, 458 

Murset, U., & Schmid, H. M. 1999, A&AS, 137, 473 

Nussbaumer, H., Vogel, M. 1987, A&A, 182, 51 

Nussbaumer, H., & Vogel, M. 1989, A&A, 213, 137 

Nussbaumer, H., Sclimid, H. M., & Vogel, M. 1989, A&A, 211, L27 

Pereira, C. B., Ortega, V. G., Monte-Lima, I. 1999, A&A, 344, 607 

Reimers, D. 1 98 1 , in Physical Processes in Red Giants, ed. I. Iben and A. Renzini 

(Dordrecht: Reidel), p. 269 
Schmid, H. M., Schild, H. 1997a, A&A, 327, 219 
Schmid, H. M. & Schild, H. 1997b, A&A, 321, 791 
Seaquist, E. R., Taylor, A. R., Button, S. 1984, ApJ, 284, 202 (STB) 
Shore, S. N., & Aufdenberg, J. P 1993, ApJ, 416, 355 
Shore, S. N., & Wahlgren, G. M. 2010, A&A, 515, A108 
Skopal, A. 2003, A&A, 401, L17 
Skopal A., 2005, A&A 440, 995 
Skopal, A. & Pribulla, T. 2006, ATel. No. 882 
Skopal, A., Pribulla, T, Budaj, J., et al. 2009, ApJ, 690, 1222 
Sokoloski, J. L., Kenyon, S. J., Espey, B. R., et al. 2006, ApJ, 636, 1002 
Tomov, N., A., Bisikalo, D. V., Tomova, M. T., & Kil'pio, E. Yu. 2010, Astron. 

Rep., 54, 628 
Vogel M., 1991, A&A, 249, 173 

Vogel, M., Nussbaumer, H., & Monier, R. 1992, A&A, 260, 156 



7 



